Non-equilibrium work fluctuations for oscillators in non-Markovian baths 
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We study work fluctuation theorems for oscillators in non-Markovian heat baths. By calculating 
the work distribution function for a harmonic oscillator with motion described by the generalized 
Langevin equation, the Jarzynski equality (JE), transient fluctuation theorem (TFT), and Crooks' 
theorem (CT) are shown to be exact. In addition to this derivation, numerical simulations of 
anharmonic oscillators indicate that the validity of these nonequilibrium theorems do not depend 
on the memory of the bath. We find that the JE and the CT are valid under many oscillator 
potentials and driving forces whereas the TFT fails when the driving force is asymmetric in time 
and the potential is asymmetric in position. 



PACS numbers: 05.70.Ln,05.40.-a 



I. INTRODUCTION 



Fluctuation theorems (FTs) which describe properties 
of the distribution of various nonequilibrium quantities, 
such as work and entropy, have been developed over the 
past decade [H, H S, 0, iH H 0] ■ Unlike most other re- 
lations in nonequilibrium statistical mechanics, remark- 
ably, the FTs are applicable to systems driven arbitrarily 
far from equilibrium. 

In particular, consider a classical finite-sized system in 
contact with a heat bath at temperature T and driven by 
some generalized force A (e.g., volume, magnetic field). 
In equilibrium, the phase space distribution of the system 
is described by a well-defined statistical mechanical en- 
semble. The free energy F can, in principle, be calculated 
and is a function of these parameters, i.e. F — F(T,X). 
At some time t, say t = 0, A is varied via a fixed path 
A(i) to a later time t = r. If this process is carried out 
reversibly, the work done W is simply the free energy 
difference AF = F(T, A T ) - F(T, A ). More generically, 
the process is irreversible and the second law gives the 
well-known inequality 



W > AF. 



(1) 



An ensemble of such processes would yield the work 
probability distribution P(W). The finite width of P(W) 
is due to two stochastic sources: (1) The system config- 
uration when the force is first initiated at t = is drawn 
from the equilibrium ensemble and (2) the path the sys- 
tem travels (in phase space) during the driving process is 
not deterministic due to the coupling with the heat bath. 
We study P(W) in the light of three closely related the- 
orems described below, the Jarzynski equality (JE) 
the transient fluctuation theorem (TFT) [1, fl 0] > and 
Crooks' theorem (CT) [1,0. 

The powerful nonequilibrium work relation due to 
Jarzynksi [3] allows one to exactly obtain equilibrium in- 
formation (the free energy difference) from measurements 
of nonequilibrium processes. The JE states 



where (3 is the inverse temperature (with Boltzmann's 
constant set to unity) and the average is over the dis- 
tribution function P(W) described above. The equality 
allows for the computation of AF even when the driving 
process is not adiabatic. This is to be compared with 
the inequality of Eq.JT]). The JE has been shown for 
Hamiltonian systems |4[ and Markovian stochastic sys- 
tems d i, Q . 

Fluctuation theorems have been found for a wide class 
of systems and various nonequilibrium quantities, includ- 
ing work, heat, and entropy production. They are also 
generally divided into steady state [3( and transient the- 
orems @, d, 0]- In this paper, we restrict our discussion 
to the transient fluctuation theorem for the mechanical 
work done. The TFT relates the ratio of probability dis- 
tributions for the production of positive work to the pro- 
duction of negative work, 



P(+W) 
P(-W) 



= e+? w . 



(3) 



Similar to the JE, the TFT has been derived in driven 
deterministic systems , and Markovian stochastic sys- 
tems |S SB laEij- 

Crooks [y] connected the JE and the TFT by using a 
relation very similar to Eq.Q, which we refer to as the 
Crooks' theorem, 



P F {+W) 
Pr{-W) 



(4) 
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>= 



-PAF 



(2) 



where Pp(+W) is the probability distribution for work 
done as described in the above scenario and Pr{—W) 
is the probability distribution for negative work done in 
a time-reversed driving process. The similarity between 
Eq.Q and Eq.(TJ| is evident and the two theorems are, 
in fact, equivalent for a large class of systems. However, 
under certain asymmetries in the potential and the driv- 
ing force, differences between the TFT and CT arise. In 
Sec. IV, we show one such scenario. 

Unlike some other nonequilibrium quantities (e.g. heat 
and entropy), the classical mechanical work is an eas- 
ily defined quantity, W = J*J" dtf(t)x(t). The original 
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formulation of the JE 4], however, relies on the gen- 
eralized work, Wj = — jl dtf(t)x(t). As discussed in 
Refs. @, (HI, these two differ by the boundary conditions 
/(t)x(t) — f(0)x(0) and care must be taken in defining 
the work in the fluctuation theorems. 

The (space and time) extensiveness of W (or Wj) 
demonstrates that second law "violations" become expo- 
nentially unlikely for thermodynamic systems; these vio- 
lations only become observably probable for microscopic 
systems. Beyond purely theoretical interests, the recent 
investigations of molecular motors and nano-mechanical 
devices demonstrate the practical importance for under- 
standing these universal nonequilibrium theorems, espe- 
cially under increasingly realistic scenarios. Technolog- 
ical accessibility of microscopic systems has opened up 
experimental study and verification of the above theo- 
rems in various systems [IJ [l2j ■ 

In this paper, we focus on the JE, TFT, and CT for sin- 
gle harmonic and anharmonic oscillators. Using second 
order Langevin dynamics of a single oscillator, the work 
fluctuation theorems have been studied in Ref . [l3[ . The 
work distribution function in harmonic polymer chains 
has been studied by one of us [14]. Both these pa- 
pers and the derivations of the FTs for stochastic sys- 
tems [E B HI S 113] utilize Mar kovian dynamics. Here, we 
study the FTs for oscillators coupled to non-Markovian 
baths. This extension is motivated by the intuition that 
there must be some memory to real heat baths and that 
all realistic models of heat baths, either classical [H| or 
quantum mechanical [lfj, are non-Markovian. Through 
analytic expressions for the harmonic oscillator and nu- 
merical simulations of anharmonic oscillators, we demon- 
strate that the memory of the bath does not affect the 
validity of the FTs. 

In the next section, we briefly describe the model and 
the generalized Langevin dynamics. Sec. Ill contains ex- 
plicit derivations of the JE, TFT, and CT in the harmonic 
limit where analytic calculations are possible. We show 
numerical results for the nonlinear oscillators in Sec. IV 
and summarize our results in the last section. 



II. THE MODEL 

We model the system as a unit mass particle in a one- 
dimensional potential with dynamics governed by the 
generalized Langevin equation 



17] 



x(t) 



dV(x) 
dx 



+ /(*)- / di' 7 (*-W) + ??(*), (5) 
o 



where V(x) is a conservative oscillator potential and/(f) 
is an externally determined time-dependent driving force. 
The rj(t) and "f(t) terms represent Gaussian noise and 
damping, respectively, and must be related through the 
fluctuation-dissipation theorem, 

<r)(t)>=0, <r 1 {t)n{t , )>=T 1 {t-t l ) : (6) 



where T is the temperature of the bath. In the white 
Gaussian noise limit, j(t) is proportional to a Dirac 5- 
function and the more familiar "Markovian" Langevin 
equation is recovered. We mention that the damping is 
even in time j(t) — j(—t), an important property to be 
used in the next section. 

We use a potential V(x) of the form 



V(x) 



k 3 x 3 



(7) 



which can represent a truncation of the Taylor expan- 
sion of some complicated potential. We restrict ourselves 
to bounded potentials, x can represent the spatial po- 
sition, an angular variable as in Ref. [l3l ]. or a gener- 
alized coordinate. From this potential and the driving 
force f(t) conjugate to x, the Hamiltonian Ti is clearly 
H(t) = x 2 /2 + V(x) - f{t)x. If the force is time- 
independent, the system would eventually reach equilib- 
rium and the free energy F can be simply calculated by 
F = —Tln(Z), where Z is the partition function. In this 
ensemble, we recall that the Jarzynski work Wj is not 
generally equal to the real mechanical work W 0, EH • In 
the next section, we derive the work distribution function 
for the harmonic oscillator, i.e. k% = k^ = 0. 



III. THE HARMONIC OSCILLATOR 

At low temperatures, many potentials can be well ap- 
proximated by the harmonic potential. Harmonic oscil- 
lators also have the practical virtue that a formal so- 
lution to Eq.(f5]) exists. We use this formal solution to 
derive the distribution functions for the Jarzynski gener- 
alized work Wj and the real mechanical work W. The 
work distribution functions are then used to verify the 
JE and the TFT. (The TFT and CT are equivalent for 
the harmonic oscillator.) The analytic expressions in this 
section are analogous to the expressions from Ref. [T3| 
with three main differences. The primary difference is 
that the noise is colored here. Second, we use the full 
second order Langevin equation instead of taking the 
strongly overdamped limit. The last difference is that 
we only study single harmonic oscillators instead of har- 
monic chains. Our results should generalize to chains, 
which are more applicable to polymer stretching experi- 
ments 12], though we do not discuss this generalization 
further. 

Equilibrium quantities are easily evaluated for the har- 
monic oscillator. Under constant driving, the free energy 
is 



F(TJ) 



Tln- 



2irT 



(8) 



Since only the free energy difference AF appears in the 
JE, only the first term on the right hand side is relevant. 
In this equilibrium ensemble, the averages for the initial 



3 



position x Q and velocity v Q and their variances are 
<x >=f/u>l, a 2 =T/u) 2 , 



< v >= 0, a 2 



T. 



(9) 



where a 2 A =< A 2 > — < A > 2 for any quantity A. 

The formal solution to Eq. ([5]) for a harmonic potential 
is (for t>0): 

x(t) = H{t)x + G{t)v + f dt'G{t-t')[f(t')+r,(t% 

Jo 

(10) 

where H (t) and G(t) are the homogeneous solutions with 
properties H(0) = G(0) = 1 and H(0) = G(0) = 0. For 
t < we define H(t) = H(-t) and G{t) = -G(-t). The 
stochastic terms are x , v OJ and r). Using the definition of 
Wj (W), we see that Wj (W) is proportional to x (x) and 
is, therefore, a linear combination of the stochastic terms. 
Thus, in the harmonic limit, Wj is Gaussian distributed 
and it is sufficient to calculate the mean < Wj > and 
variance a 2 . Obviously, W is Gaussian as well and its 
distribution function is 



P(W) 



1 



(W-<W» A 



V27TCT 2 



(11) 



For such Gaussian processes, the TFT is satisfied if < 
W >= f3a 2 / 2 and the JE is satisfied if < Wj >= AF + 
f3o 2 /2. 

In order to compare the means and variances of the 
work distributions, we first derive identities for the 
Green's functions H(t) and G(t) in Eq. (fTOf . The Laplace 
transform of Eq. ([5]) is 

s 2 x(s) — sx — v a + j(s)[sx(s) — x a ] + lo1x{s) = 0, (12) 

where x(s) = dtx(t)e~ st is the standard definition of 
the Laplace transform. The Green's functions, H(t) and 
G(t), must also satisfy equations analogous to Eq. (fT^|) 
and simplify due to their initial conditions. Solving these 
algebraic expressions gives 



H(s) 



7 + s 



S 2 + + LU 2 ' 



G(s) = -= = ? . 



(13) 



Some manipulations and the inverse transform reveal 
identities between the two Green's functions in Laplace 
and real space, 

sH(s) = l-^G(s), 
sG(s) = H(s)-i(s)G(s), 
H{t) = -u 2 G(t), 

G{t) = H{t)- [ dt'i(t-t')G(t'). (14) 
Jo 

For white noise, where H (t) and G(t) are well-known, we 
confirm that Eqs. fT^)) are correct. 



We plug Eq. (fT0"f into the definition of the Jarzynski 
work (for a time r) , 

Wj = - [ dtf(t)[H(t)x + G(t)v ] 



dt / dt'f(t)G(t-t')[f(t') + V (t')]. (15) 
o Jo 



With the use of the equilibrium averages, we find the 
mean of the Jarzynski work, 



<Wj> 



dtf(t)H(t)f(0)/u; 2 o 



dt / dt'f(t)G(t-t')f(t'). (16) 
o Jo 



We can re-express Eq. (|16l) for later comparison with a 2 
by integrating by parts and using the identity for H(t) 
in Eq.lJTJ" 



< Wj >= AF + / dt dt'f(t)H(t - t')f(t')/uj 2 . 
Jo Jo 

(17) 

We find the variance of Wj by using the oscillator equi- 
librium averages and the fluctuation-dissipation theorem, 



P°j = 



dtf(t)H(t) 



1 



dtf(t)G(t) 



dh / dt[ / dt 2 / d^7(<i - 1' 2 ) 
Jo Jo Jo 

x /(ti)G(*i-ti)/(t a )(?(t 2 -^). (18) 

In order to simplify this expression to compare with 
Eq. (|17p , we must reduce the quadruple integral to a more 
manageable double integral. We define 

i(h,t 2 )= r dt[ f 2 dt^h ~ - t 2 h(t[ - 1 2 ). 



(19) 

The integrals of I(ti,t 2 ) can be evaluated by first 
doing a double Laplace transform, I(si,s 2 ) — 
]™dti J °° dt 2 e- Sltl e- S2t2 I(t 1 ,t 2 ). This double trans- 
form can be done by using the even symmetry of ^(t). 
We separate and define the two symmetric parts of "f(t) 
using the step function, 

j(t) - 7+ (t) + 7 _(t) = j(t)Q(t) + j(-t)Q(-t). (20) 

Incidentally, the Laplace transforms of the two separate 
parts are equal to the Laplace transform of 7(t), 7+(s) = 
7_(s) = 7(5). We use this property, the convolution 
theorem, and Eqs. fT^)) to find the double transform, 

- H( Sl )+H( S2 ) H( Sl )H(s 2 ) ~ ~ 
I{si,s 2 ) = j- ■ 5 G(si)G(s 2 ). 



U 2 {S! + S 2 ) 



(21) 
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The inverse transform can easily be done, giving 
Wuh) = ggi^, -G( tl )G(t 2 ). (22) 

U! LU 

o o 

Finally, we plug this expression back into Eq. I|18p , 



(3<Tj = dt dt'f(t)H(t~t')f(t')/oj 2 
Jo Jo 

= 2<Wj> -2AF. (23) 

The last equality proves the Jarzynski equality for har- 
monic oscillators, even when r)(t) is not (5-correlated. As 
in Markovian stochastic derivations of the TFT @, 0, [H[ , 
the generalized work Wj does not satisfy the TFT, how- 
ever Wdiss — Wj — AF does. 

A simpler derivation of the JE follows if we assume 
time-translation invariance of various correlation func- 
tions. The definition of a 2 contains the auto-correlation 
of Ax(t) = x(t)- < x(t) >, 

aj = f dt f dt'f(t)f(t') < Ax(t)Ax(t') > . (24) 
Jo Jo 

We evaluate the correlation in the brackets by using time- 
translation invariance and the formal solution of Ax(t) = 
H(t)Ax + G{t)v a + J* dt'G(t - t')r)(t'), 



< Ax(t)Ax(0) >= H{t)/{(3u 2 ) 



(25) 



Plugging this back into Eq. ([24|) . we recover the result 
from the previous derivation Eq. ([23)) and the JE is easily 
seen. 

A similar analysis is done for the real mechanical work, 



W = 



dtf(t)[H(t)x + G(t)v 



+ dt dt'G{t-t')[f(t')+>n{t')]. (26) 
Jo Jo 

The average work is obtained using the equilibrium av- 
erages, 

<W> = f dtf(t)H(t)f(0)/to 2 
Jo 

+ f dt f dt'f(t)G(t-t')f(t'). (27) 
Jo Jo 

The Laplace transform manipulations and the symmetric 
property of j(t) can be used to find an expression for the 
variance of the work a 1 . However, we only show the 
simpler derivation, analogous to Eq. (|2"4"l) , 



dt \ dt'f{t)f{t') < Ax(t)Ax(t') >, (28) 




sinusoidal 
linear 

symmetric sawtooth 
asymmetric sawtooth 



where Ax(t) = x{t)— < x(t) > . The quantity in brack- 
ets can be calculated by using the formal solution for 
the velocity, i.e. the time derivative of Eq. (TlTjj) . Timc- 



FIG. 1: (Color online) The inset shows the work probabil- 
ity distribution P(W) for a harmonic oscillator under a driv- 
ing force f(t) = sin(nt/T) and a bath with exponentially 
correlated noise. Simulation details are given in Sec. IV. 
The distribution is a Gaussian with a positive mean, but 
with significant amplitudes for negative W. Work probability 
distributions for other driving forces are qualitatively simi- 
lar, as is P(Wdiss)- The main figure shows a log-lin plot of 
P(W) / P(—W) versus W. The different points are for a va- 
riety of driving forces and their excellent agreement with the 
exponential confirms the validity of the TFT. 



translation invariance is again assumed and the velocity 
auto-correlation is easily calculated 

< Ax(t)Ax(0) >= G(t)/f3. (29) 

We use this result in Eq. (f2"5)l and find 

(5a 2 = f dt f dt' f{t)G(t-t')f(t') 
Jo Jo 

= 2 < W > . (30) 

The last equality is from a comparison with Eq. (|27p and 
is valid when /(0) = 0. Under this condition, we thus 
prove the TFT for the probability distribution for the 
real mechanical work. 

Eq. (j2"3")) and Eq. (f3"0|) are the main results of this sec- 
tion. Simulations of a driven harmonic oscillator in a non- 
Mar kovian bath confirm these derivations; for all driving 
forces and bath conditions simulated, the JE and TFT 
are true for the harmonic oscillator. Figure [T] shows the 
results of these simulations. Details of the numerics and 
simulation results for anharmonic oscillators are given in 
the next section. 



IV. ANHARMONIC OSCILLATORS 

For anharmonic oscillators, finding simple expressions 
for P(W) and P(Wj) is, in general, not possible; W and 
Wj are not linear combinations of the Gaussian stochas- 
tic quantities, therefore their distributions are not Gaus- 
sian. In this section, we numerically measure the work 
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distribution functions for oscillators with the potential of 
Eq.flTJ) when k 3 , k A > 0. 

Random uncorrelated Gaussian numbers are com- 
monly numerically generated and used in standard in- 
tegration algorithms for Langevin equations of motion. 
An efficient Verlet-like integrator for second order white 
noise Langevin equations is given in Ref. 18] . We mod- 
ify this algorithm in order to simulate Langevin systems 
with exponentially correlated noise. The damping kernel 
has the form 



7(*) 



-r|*| 



(31) 



Exponentially correlated noise can be effectively intro- 
duced by using white noise in the equation of motion for 
an auxiliary variable z, in addition to the position x and 
velocity v. The coupled differential equations are 



dx/dt = v, 

dv/dt = -dV(x)/dx + f(t) + z, 
dz/dt = — Tz — v + £. 



(32) 



C represents Gaussian white noise which is easily gener- 
ated numerically, 

<C(*)>=0, <((t)((t') >= 2T5(t-t'). (33) 



It can be shown that the equations of motion Eas. (|32)l are 
equivalent to Eq.([5]) with exponentially correlated noise 
Eq . (f3"Tj) . In all simulations, T = T = ui^ = 1, though we 
have verified that our results are qualitatively the same 
with different temperatures and damping constants. We 
use a step size h with h = 0.01 in all simulations shown. 
We have checked that our results do not change with 
smaller step sizes. 

Work distributions are approximated by histograms of 
O(10 6 ) measurements of the work done (W and Wj) over 
a time r. Between each measurement, we allow the os- 
cillator to equilibrate by integrating Eqs. (f3"2")l for O(10 5 ) 
steps with fit) = 0. No differences in the distributions 
are seen for different histogram bin sizes and longer equi- 
libration times. 

We use a sawtooth driving force of the form, 



f{t) = t/t , 0<t<t o 
= i-t + T)/(r-t ), 



t <t<r. (34) 



Under this driving force, the equilibrium configurations 
at t = and t = r are identical, i.e. AF = and 
W = Wj. In the r — » limit, the work distribution is 
trivially peaked at W = because there is no time for 
work to be done. P(W) is similarly peaked at zero in the 
t — > oo limit because the system is driven adiabatically 
from an equilibrium state back to the same equilibrium 
state. Our numerical simulations are in accord with these 
physical limits. The figures show simulation results with 
r = 10, which is intermediate between these two lim- 
its. (No qualitative differences in terms of the JE, TFT, 
and CT exist with different r.) Lastly, /(0) = 0, which 




FIG. 2: (Color online) The inset shows the work probability 
distributions for oscillators with quartic anharmonicities un- 
der sawtooth driving with t /r — 0.25, 0.5, and 0.75. These 
distributions have positive means, but also substantial weight 
on the negative region of W. Po.25t(W) = Po.75r(W), indi- 
cating that the asymmetry in the driving force does not alter 
the work distribution functions, thus the ratios in the defini- 
tion of the TFT and the CT, Eq.lpH)) and Eq.fl36j, must also 
be equal. The main figure confirms this equality and shows 
that the TFT and the CT are valid for the quartic oscillator. 



is a necessary condition for the validity of the TFT for 
the harmonic oscillator detailed in the last section and 
Figure [T] 

By changing t a , we can alter the symmetry of fit); 
the force is symmetric in time only when t Q = 0.5r. We 
implement this symmetric driving and two asymmetric 
sawtooth forces with t a — 0.25r, 0.75r in our simulations. 
As mentioned in the introduction, there is a subtle dif- 
ference in the TFT and the CT with the latter using a 
time-reversed process in the denominator of the ratio of 
probabilities. For the symmetric sawtooth force, forward 
and reverse driving must be equivalent and no differences 
exist between the TFT and the CT. Our asymmetric saw- 
tooth forces are complementary in that the time-reverse 
driving of one force is equivalent to the time-forward driv- 
ing of the other. In other words, in terms of this sawtooth 
force, the TFT is 



PtS+w) w 

Ptoi-W) 



whereas the CT is 



PtJ+W) 

Pr-t i-W) 



(35) 



(36) 



(P S {W) is the work distribution corresponding to a saw- 
tooth potential with a break at time t = s.) 

We first simulate an anharmonic spring with a unit 
quartic anharmonicity (i.e. k 3 = 0,^4 = 1). Figure [2] 
displays the results of these simulations. Aside from not 
being Gaussian, the probability distributions have similar 
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For our example the validity of the JE follows from the 
CT. Independently, we have also measured the average 



exponential of the work and find < e 



-pw 



>= 1 within 



o 
w 



TFT t o =0.25x 
TFT t o =0.75x 
10 ' TFT.CT t o -0.5x 
CT t„=0.25x 



0.6 
W 



FIG. 3: (Color online) The top panel shows the work prob- 
ability distributions for particles in an asymmetric poten- 
tial under sawtooth driving with t jr = 0.25,0.5, and 0.75. 
Due to the lack of symmetry in the potential, Po.25t{W) 
Po.75t{W). The bottom panel shows the ratios from Eq. (|35p 
and Eg . (I36fl versus W on a log-lin plot. The data sets that 
agree with the exponential indicate that the CT is valid even 
for asymmetric potentials. The other sets indicate that the 
TFT fails for these potentials under general driving forces. 



properties as the work distribution of the harmonic oscil- 
lator: the means are all positive and there is a very sub- 
stantial probability of measuring negative work. Further- 
more, the two complementary asymmetric driving forces 
give the same work probability functions. From Eq. (|35j) 
and Eq. ([36|) . it is clear that, in this case, the TFT and 
the CT are equivalent. The main plot in Fig. [2] displays 
this equality as well as the validity of the theorems. 

When a cubic nonlinearity is included (£3 = fct = 1), 
the even symmetry of the potential is broken. Figure [3] 
displays results for these non-symmetric oscillators. The 
top plot is analogous to the inset of Fig. [2] and shows 
the probability distributions for the three different saw- 
tooth driving forces. From this figure, it is clear that all 
three distributions are different and the ratio in the TFT 
Eq. (f3"5)) is not equivalent to the ratio in the CT Eg . (f3T>)) . 

The difference between the TFT and the CT is clearly 
shown by the bottom panel of Fig. [3l In fact, the CT 
is valid under all sawtooth driving forces, whereas the 
TFT is only valid in the symmetric driving force. For 
oscillator potentials that are not symmetric in x, i.e. 
V(x) ^ V(—x), and the driving force is asymmetric in 
time about r/2, the TFT dramatically fails. 



~ 4 decimal places for both symmetric and asymmetric 
oscillators. Because AF = 0, this indicates the Jarzynski 
equality Eq.@ is most likely valid for a large class of 
oscillator potentials, even with exponentially correlated 
noise. 

When a sinusoidal force of the form f(t) = sin{mct/r), 
where n is an odd integer, is used, the simulation re- 
sults are qualitatively the same as with the symmetric 
sawtooth force. Obviously, this is because the period of 
the sinusoidal force gives a symmetric (about r/2) driv- 
ing process. We also simulate linear driving which takes 
the system to a different equilibrium configuration in the 
adiabatic limit. Though AF 7^ with linear driving, we 
find that the JE and the CT are still valid, but the TFT 
fails. This last result is not surprising because the lin- 
ear driving force has many qualitative similarities with 
asymmetric sawtooth driving forces. 

The above results are unchanged when colored noise is 
replaced by white noise. The white noise results are to 
be expected due to the numerous analytical derivations 
of these nonequilibrium theorems [1, H, 0, H, H, EH for 
Markovian baths. Our simulation results shown in this 
section indicate that the CT and the JE are also valid 
for non-Markovian baths. 



V. SUMMARY 

We have derived non-equilibrium work fluctuation the- 
orems for the classical harmonic oscillator connected to 
a generalized Langevin bath by studying the work dis- 
tribution functions of the real mechanical work and the 
Jarzynski generalized work. Both W and Wj are Gaus- 
sian variables. We derive the TFT for the real work and 
the JE by using exact relations between < W > and a 2 
and < Wj > and Oj. To our knowledge, these are the 
first rigorous derivations of the nonequilbrium work fluc- 
tuation theorems for a system described by an arbitrary 
damping kernel of the generalized Langevin equation. 

We also numerically measure the work distribution 
functions for anharmonic oscillators in Langevin baths 
with exponentially correlated noise. Our numerical sim- 
ulations indicate that the JE extends to particles in a 
variety of oscillator potentials and driving forces. They 
also show that the symmetries of the potential and the 
driving force account for differences between the TFT 
and the CT. Crooks' theorem is valid for all potentials 
and driving forces simulated, whereas the TFT is not 
valid when the oscillator potential is not even in x and 
the driving force is not symmetric in t (about r/2). The 
strong agreement between our simulation results and the 
nonequilibrium theorems (CT and JE) provides motiva- 
tion for a derivation of the theorems for anharmonic os- 
cillators in non-Markovian baths. 

These results are important due to the existence of 
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noise correlations in any real heat bath. This has obvi- 
ous consequences for experimental tests of the Jarzynski 
equality and the fluctuation theorems. Experimentally, 
the JE is a powerful tool to measure equilibrium free 
energy differences due to the fact that any real driving 



is done irreversibly. Interesting and open questions still 
exist on the ramifications of the noncquilibrium fluctua- 
tion theorems and the second law "violating" events for 
molecular engines and microscopic thermodynamics. 
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